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ABSTRACT 

We investigate the spatial clustering of dark matter halos, collapsing from 1 — 4cr 
fluctuations, in the redshift range 0 — 5 using N-body simulations. The halo bias of 
high redshift halos {z ^ 2) is found to be strongly non-linear and scale-dependent on 
quasi-linear scales that are larger than their virial radii (0.5 — 10 Mpc/h). However, 
at lower redshifts, the scale-dependence of non-linear bias is weaker and and is of 
the order of a few percent on quasi-linear scales at 0 ~ 0. We find that the redshift 
evolution of the scale dependent bias of dark matter halos can be expressed as a 
function of four physical parameters: the peak height of halos, the non-linear matter 
correlation function at the scale of interest, an effective power law index of the rms 
linear density fluctuations and the matter density of the universe at the given redshift. 
This suggests that the scale-dependence of halo bias is not a universal function of 
the dark matter power spectrum, which is commonly assumed. We provide a fitting 
function for the scale dependent halo bias as a function of these four parameters. Our 
fit reproduces the simulation results to an accuracy of better than 4% over the redshift 
range 0 ^ z ^ 5. We also extend our model by expressing the non-linear bias as a 
function of the linear matter correlation function. It is important to incorporate our 
results into the clustering models of dark matter halos at any redshift, including those 
hosting early generations of stars and galaxies before reionization. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - galax¬ 
ies: statistics - galaxy: haloes 


1 INTRODUCTION 

The spatial distribution of luminous galaxies is a valuable 
resource for probing cosmology and the physics of galaxy for¬ 
mation. The clustering of the galaxy distribution is shaped 
by the clustering of the dark matter halos which host them. 
The clustering of dark matter halos can be quantified us- 


the clustering of halos using scale independent, linear bias is 
expected to be valid on scales larger than the virial radii of 
halos where dark matter halo substructure is not important. 

However, the simple picture of a scale-independent halo 
bias has been shown to be inaccurate and various non-linear 
and non-local processes result in some degree of scale depen- 


in a non-linear and deterministic fashion (Fry & Gaztanaga 


ing the halo bias which describes how dark matter halos dence 

1 Matsubara||1999 |Cole et al. 

2005| Seo k, Eisenstein| 

trace the dark matter 1 Kaiser 1984 Bardeen et al. 1986 

20051 

Angulo et al. 

20051 Huff et al. 

20071 ISmith et al.|20071 

Bond et al. 19911. Conventional models assume that the halo 

Angu 

o et al.|20081 

McDonald & Roy|20091 Desjacques et al.| 

bias is related to the underlying dark matter density field 

2010F 

Musso et al.|2013||Paranjape et al.|2013|. Incorporat- 


1993). Mo & White (1996) showed that, on large scales, the 


halo bias can be approximated as a scale independent func¬ 
tion of the mass of the halos. In particular, they showed 
that the clustering of dark matter halos is proportional to 
that of the dark matter with the constant of proportional¬ 
ity being called the linear halo bias. The approximation for 
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ing such a scale dependence of halo bias into theoretical 
models could be crucial for interpreting the clustering of 
galaxies. While several studies have focussed on the scale 
dependence of the bias on very large scales, its scale depen¬ 
dence on scales larger than the typical virial radii of dark 
matter halos is equally interesting. These scales, correspond¬ 
ing to comoving length scales of 0.1 to a few megaparsecs, 
are smaller than scales where the matter distribution is still 
linear and therefore are referred to as quasi-linear scales. The 
scale-dependence of halo bias on these scales arises mainly 
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due to the non-linear growth of matter fluctuations (Smith 


et al. 20071 and is difficult to estimate using perturbative 


approaches because of the non-linearity of matter density 
held ( Reed et al.|2009 |. 

There have been studies in the literature of deviations 
from the linear bias approximation on quasi-linear scales us¬ 


ing analytic techniques (Scannapieco &: Barkana|2002 Iliev 


et al.|2003 Scannapieco fc Thacker|2005 I as well as N-body 


simulations (|Hamana et al.||2001| |Diaferio et al.||2003| |Cen| 

et al. 1120041 |Tinker et al. 2005[ |Gao et al. |2005b|' Angulo 


et al.|2008 Reed et al.|2009 van den Bosch et al. 12013 1. In 

particular, these studies focussed on the clustering of dark 


2.1 The linear bias model for halo clustering 

In the linear bias approximation, the cross-correlation be¬ 
tween halos of mass M' and M" is given by 

6.h(r|M', M", z) = h{M\ z)b{M", z)^^^{r, z), (1) 

where z) is the non-linear two point correlation func¬ 

tion of matter density contrast at redshift 2 and b(M, z) 
is the scale independent linear bias of halos of mass M at 
this redshift. Eq. 0 is valid on large scales, where density 
perturbations grow linearly with redshift (Cooray & Sheth 


20021. The two-point matter correlation function is obtained 


matter halos either in the local universe (eg. Tinker et al. 


2005 1 or at very high redshifts before the reionization of the 
intergalactic medium ( |Reed et al.||200'^ . In general these 
studies showed that the halo bias is non-linear and scale- 
dependent on quasi-linear scales, but the scale dependence 
weakens on large scales. Specifically, |Reed et al.| ( |2009| ) find 
a strong scale dependence of halo bias on quasi-linear scales 


by Fourier transforming the non-linear matter power spec¬ 
trum, P(k,z) ( Smith et al.|20d3 l 


OO 

^) = j ^ 


sin(A:r) 

kr 


( 2 ) 


It is well known that the scale independent halo bias can 


for rare dark matter halos at high redshift, with the scale 

be expressed as a 

function of the ‘peak height’, v{M,z) = 

dependence increasing with the rarity of the halo. 

Sc/cr{M, z), of dark matter halos ( 

Mo & White|1996| 

Sheth 

The motivation of this paper is to study the cluster- 

^ Tormen||1999[ ISheth et al.||2001| |Cooray & Sheth 

|2002| 

ing of dark matter halos with a specific focus on the scale 

Tinker et al.|2010 

1. The peak height is a measure of the rar- 

dependence of halo bias on quasi-linear scales. In particu- 

ity of halos (jSheth et al. 12001 1 with rarer halos having larger 


knowledge, no such previous studies have been carried out. 
This will help to gauge the amplitude, scale dependence and 
evolution of the bias of dark matter halos for 0 sC a ^ 5 and 
hence bridge the gap between other studies which focus on 
the epochs before reionization. We will address this issue us¬ 
ing N-body simulations to measure the dark matter and halo 
correlation functions in the real space. These measurements 
will be used to calibrate the nature and evolution of the 
non-linear halo bias in the redshift range 0 — 5 over a range 
of length scales. In particular, we find that the bias of dark 
matter halos is non-linear and scale dependent on quasi- 
linear scales. Furthermore, it is not possible to express this 
scale dependence in terms of the usual parameterizations 
and therefore one has to invoke additional parameters. 

The organization of this paper is as follows. In Section 2, 
we compare the halo bias of rare high redshift dark matter 
halos computed from analytic models and simulations. In 
Section 3, using simulations, we probe the scale dependence 
and redshift evolution of the non-linear bias of rare halos 
in the redshift range 0 to 5 and obtain a htting function to 
describe these effects. We conclude with a brief discussion 
of our results and their implications in the final section. 


^{M, z). Here, Sc = 1.686 is the critical density for halo col¬ 
lapse and cr(M, z) is the rms linear density fluctuation on a 
mass scale M 


a\M,z) = a^{R,z) = j ^k^P'^’^{k, z)W\k, R), 
0 


(3) 


where R is the comoving radius of a sphere containing mass 
M, W{k,R) is the Fourier tranformation of the top hat win¬ 
dow function and P^^"{k, z) is the linear matter power spec¬ 
trum. 

In particular, for linear halo bias, we use the fitting 


function of Tinker et al. (20101 which was calibrated against 


N-body simulations and is given by, 
b[M,z) = b[v{M,z)) = l- A 


+SS 


+ bJ' + Cu\ (4) 


I Tinker et al.| {[20To| estimate the free parameters of Eq. |4| 
to be T = 1.0, a = 0.132,5 = 0.183,6 = 1.5, C = 0.265 
and c = 2.4. The halo bias given by Eq. (|^ increases with 
increasing u{M,z). 

We assume that dark matter halos which can host galax¬ 
ies have a spherical over-density A = 200 times the average 
density of univers^ Then, the virial radius r 2 oo of a halo of 
mass M is M = (4/3)7rrfooPcA. Under the above assump¬ 
tions the halo correlation function is 


2 CLUSTERING OF RARE DARK MATTER 
HALOS AT HIGH-Z 

In this section, we investigate whether the linear bias model 
for halo clustering gives a good description of the clustering 
of high-z dark matter halos on quasi-linear scales. For this 
we first describe the linear bias model for halo clustering in 
Section [2.1 1 In Section [2.2| we introduce the N-body simula¬ 
tions used in our study. The clustering of dark matter halos 
estimated from these simulations is then compared with the 
clustering prediction using the linear bias model. 


1 + ^hh(r|M', M", z)=[l + 6(M', z)b{M", «)^mm(r)] 

e[r-r,^in(M',M")], (5) 

where the function 0[r — rniin{M', M")\ incorporates halo 
exclusion to ensure that , M", z) = —1 for rmin = 

max[r 2 oo(AF'), r 2 oo{M")]. 


^ [Tinker et al.| | |2010| calibrate their fitting function for the large 
scale bias as a function of A. Here, the quoted parameter values 
are for halos with A = 200. 
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Figure 1. Upper panels : The two point correlation functions of dark matter halos in the mass range ~ IQl'^ — 5 X at various 

redshifts as labelled. The points (triangles and circles) are measured from N-body simulations and the curves (solid and dotted) are 
analytic predictions using the linear bias approximation with the same cosmological parameters as used in the simulations. The results 
at 2 : = 2 are from the MXXL simulation and those at other redshifts are from the MS-W7 simulation. Bottom panels: the ratio of the 
correlation functions measured from simulations to those computed analytically. 


The two point correlation function of dark matter halos 
in a mass bin M' ^ M ^ M" is given by 


l + ^hh(r|[M^M^,^) = 


A1 IVl 

U) / "■ / 


dM2 


, M" 

M' M' 

n{M')n{M") [1 + Chh(r|Mi, M 2 , 2 )], (6) 

where n{[M', M"], z) = jff, dMn{M, z) is the total num¬ 
ber density of halos in the mass bin [M', M”]. For the halo 


mass function, n(M, 2 ), we use the fitting function of Jenkins 
|et al.| ( |2001[ ) which is in excellent agreement with the mass 
functions obtained from the simulations used in our study. 
Eq. § is a reasonable approximation for the usual 2-halo 
term for halo clustering on scales larger than the virial radii 
of dark matter halos ( Cooray fc Sheth||2002 K 

The average bias of halos with mass between M' and 
M”, on scales bigger than their virial radii can be written 
as being scale-independent and is given by 


M" 


6([M',M"],2) = 


n([M', M"], z] 


j dMb{M,z)n{M,z). 


Using this result for the bias in Eq. Q, we get 


( 7 ) 


Chh(r|[M',M"], 2 ) = b^{[M',M"],z)^^^{r,z). (8) 

In what follows, we will compute the halo correlation func¬ 
tions for dark matter halos in mass bins and compare with 
those measured from N-body simulations. 


2.2 N -body simulations 

Our study mainly uses two cosmological dark matter N-body 
simulations, the MS-W7 simulation ( |Guo et al.||2013] |Pike| 
et al.||2014 1 and the Millennium-XXL or MXXL simulation 
(Angulo et al. 20121. The MS-W7 simulation uses a cubic 


computational box of comoving length 500 h~ Mpc with 
2160® particles of mass 8.61 x 10®M©. This is used to probe 
the clustering of halos at 2 = 3, 4 and 5. It adopts a fiat 
ACDM background cosmology, which is in agreement with 
the WMAP7 results (Larson et al. 20111, with h = 0.704, 


Qb = 0.0455, Qc = 0.2265, = 0.0, erg = 0.81 and m = 

0.967. 

The MXXL extends the previous Millennium and 


Millennium-II simulations (Springel et al. 2005 


Boylan- 

Kolchin et al.|2009 1 and follows the evolution of 6720® dark 
matter particles inside a cubic box of length 3000 h~^ Mpc. 
The particle mass is 8.46 x 10®M©. This simulation adopts a 
ACDM cosmology with the same cosmological parameters as 
the previous Millennium simulations. Accordingly, h = 0.73, 
fib = 0.045, flc = 0.205, fli, = 0.0, erg = 0.9 and ris = 1.0. 
The MXXL halos are used to investigate the clustering at 
2 = 0,1, 2 and 3. We also compare the results obtained using 
MXXL simulation at 2 = 3 with those obtained using using 
the Millennium simulation at the same redshift, which has 
a box of 500 h~^ Mpc. 

In both simulations, groups of more than 20 particles 
are identihed as dark matter halos using a friends-of-friends 
algorithm (FOF(0.2)) with linking parameter equal to 0.2 of 
the mean particle separation (Davis et al. 19851. The halo 


mass functions from these simulations are well described by 
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4 C. Jose et al. 


the fitting function of Jenkins et al. (20011 over a wide range 
of halo masses and to an accuracy better than 10 %. 

The two point correlation functions of dark matter ha¬ 
los and dark matter particles from the simulations are com¬ 
puted by counting the number of pairs as a function of the 
separation, r, relative to that of a random distribution and 
is given by 


C (r) 


N^{r) 

N?an{r) 


(9) 


where Np{r) is the total number of pairs in the simulation 
separated by a distance r to r -|- Jr and is the total 

number of pairs. 

As mentioned above, in this paper, we focus on the clus¬ 
tering of rare dark matter halos on quasi-linear scales in the 
redshift range 0 — 5. Therefore, we consider only those ha¬ 
los with a peak height u{M,z) > 1. At 2 : ~ 0, the typical 
masses of these halos range between 10^^ — W^^Mq and 
therefore they correspond to poor galaxy groups and clus¬ 
ters. On the other hand, for z ^ 2, the masses of rare halos 
range from 10^° — IO^^Mq. As we see later, the scale depen¬ 
dence of the halo bias due to non-linear clustering is much 
more signihcant at higher redshifts (2 = 2 — 5) than it is at 
lower redshifts. Therefore, we first address the issue of the 
non-linear clustering of high redshift dark matter halos on 
quasi-linear scales and then its evolution in the low redshift 
universe. 



Figure 2. The halo bias b{r, M, z) = com¬ 

puted at 2 = 3 using Millennium and MXXL halos in z) 
bins, as indicated by the label. The thick (blue) and thin (black) 
lines respectively correspond to halos with 2.41 v ^ 2.28 and 
3.07 ^ V ^ 2.72. The thick (blue) and thin (black) vertical lines 
indicate scales corresponding to twice the virial radius of the most 
massive halo in each sample. 


2.3 Comparing simulations and linear bias models 

We hrst show that the clustering strength of high-z, rare 
dark matter halos on quasi-linear scales differs significantly 
from the predictions of the linear bias model by compar¬ 
ing with the spatial correlation functions estimated from N- 
body simulations. In the top panels of Fig. the halo cor¬ 
relation functions estimated from simulations r, 2 )) 

are shown at 2 = 2, 3, 4 and 5 for halos in the mass range 
9 X 10^° — 10^^ (black circles) and 2 x 10^^ — 4 x 10^^ 

h~^M q (red triangles). We note that, these halos respec¬ 
tively host typical Lyman-a emitters (LAEs) and Lyman- 
break galaxies (LBGs) in the same redshift range 
2013b|a| and are rare halos, collapsing from 2 — 
ations {y{M,z) ~ 2 — 3). On small scales the correlation 
functions drop to —1 due to halo exclusion. These scales 
correspond to the typical virial radius of halos in the given 
mass bin. 

Also shown in the top panels of Fig. [^are the correla¬ 
tion functions, ^hhiM, r, 2 ), predicted by Eq. ([^ for the same 
cosmological parameters as used in the simulations. The cor¬ 
relation functions are computed for halos in the same mass 
bins used to estimate r, 2 ). Fig. [^clearly shows that 

2 ) and ^hh(r, 2 ) agree well with each other on large 
scales (r > 10 — 15 h~^ Mpc). However, on quasi-linear scales 
(r ~ 0.5 — 10 h~^ Mpc), determined from the sim¬ 

ulations shows an excess compared to ^hh(r, 2 ) computed 
using Eq. §. 

To understand the degree of this deviation more clearly, 
we have plotted in the bottom panels of Fig. [^the ratio of 
the dark matter halo correlation functions measured from 
simulations to that computed from the linear bias model 
(i.e ■*)/^hh(r, 2 )) for each mass bin. It is clear from 


Jose et al. 
3 (t fluctu- 


the figure that, on quasi-linear scales, the predictions of the 
scale-independent bias model are insufficient to explain the 
halo correlation functions measured directly from the simu¬ 
lations. For example, the massive halos at the highest red¬ 
shift (2 X 10^^ ^ M/Mq 4 X 10^^ at 2 = 5) show clustering 
in the simulations that is sometimes larger by a factor as 
large as ~ 20 at 1 ^ r ^ 2 Mpc/h, compared to the linear 
bias model predictions. On the other hand, at lower red¬ 
shifts and for less massive halos (9x 10^° ^ M/Mq ^ 10^^ at 
2 = 3), the clustering excess is only a factor of 2—3 at r ~ 0.5 
Mpc/h. Furthermore, the deviation between Chi“(r, 2 ) and 
^hh(f', z) increases with the redshift and mass of dark matter 
halos. Overall we conclude that the halo bias of high redshift 
dark matter halos is strongly scale dependent on quasi-linear 
scales and the scale dependence increases with the rarity of 
the halos. 

Earlier studies focused on the non-linear bias of halos at 
the present epoch ([Hamana et al.|2001[ |Diaferio et al.|2003| 


Tinker et al.|2005 I or at redshifts before reionization ( |Reed 

et al. 20091. The scale dependence of the non-linear bias 
in the fitting function s provided by|Hamana et al.| ( |2001[ ); 
Diaferio et al. (20031; Tinker et al. I (20051 is too weak to 


explain the clustering of halos at the redshifts, masses and 
scales of interest here. The fitting function of |Reed et al.| 
(20091 has a stronger scale-dependence and describes the 


non-linear clustering of high redshift MS-W7 halos correctly. 
However, as we see later, their results are not consistent with 
the bias measured from the MXXL simulation and also at 
lower redshifts (2 = 0 — 2). Therefore, non-linear clustering 
of rare halos on quasi-linear scales has not been satisfactorily 
addressed and thus warrants further investigation. 
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The clustering of dark matter halos: scale-dependent bias on quasi-linear scales 5 



Figure 3. The halo bias b{r,M^z) = yj plotted as a function of the redshift range 2 — 5. Each panel 

shows the b{r, M, z) of halos in ^{M, z) bins, as indicated by the label. The results obtained from the MXXL simulation are shown by 
the solid black lines {z = 2) and red circles {z = 3) whereas other curves corresponds to the results from the MS-W7 simulation. 


3 THE SCALE-DEPENDENT, NON-LINEAR 
HALO BIAS 

3.1 The measured bias 

We have shown in the previous section that, on quasi-linear 
scales, high-z dark matter halos collapsing from 2 — 3 ct fluc¬ 
tuations cluster more strongly than the predictions of the 
linear bias model. Therefore, to understand the clustering 
of these rare halos, one has to invoke a scale-dependent, 
non-linear bias. For this, we first define a non-linear, scale 
dependent halo bias of dark matter halos at any redshift as 


( 10 ) 

Here, is fh® non-linear dark matter correlation 

function computed directly from the simulations using 
Eq. The function b(r, M, z) is thus expected to be inde¬ 
pendent of r on large scales. 

We note that there are alternative definitions of the 


i Scannapieco fc Barkana]|2002 Reed et al.|2009 


b„i{r, M, z) = 


ehr(r,g) 


halo bias in Fourier space and also as the ratio of the halo- 
matter cross correlation function to the matter correlation 
function (eg. Tinker et al.|2010 Manera & Gaztanaga 20111. 
The choice of a particular definition of the halo bias can in 
principle introduce a scale-dependence in the bias measured 


from the simulations (Baumann et al. 20131. This scale- 


dependence introduces a few percent difference between the 


measured bias that is defined by Eq. (101 and those de¬ 
fined by other bias definitions ( Smith fc Marian||2011 Pol¬ 
lack et al. 2012 b This effect is much weaker compared to 


the strong scale-dependence of the halo bias presented in 
this work and hence will be neglected. We also note that 
bni{r, M, z), defined by Eq. (101, encapsulates the scale- 
dependence of the halo bias due to the non-linear higher 
order correlations of the matter distribution and also due to 
the scale-dependence of the linear halo bias. Thus, one can 
directly use (r, M, z) to compute the ^hh at any given 
scale in real space using Eq. |5f with minimum ambiguities, 
which is our primary goal. 

In the previous section, we estimated the clustering of 
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(^min ? t^max) 



MXXL halos 

MS-W7 halos 

Millennium halos 

z= 

=0 

Z=1 

z=2 

z=3 

z=S 

Z = A 

z=5 

z=S 

(3.37,4.21) 

3.77 

1.9 X 

IO13 

3.0 X loll 

6.1 X IO13 

1.5 X 1013 

8.0 X 10l3 

2.7 X 10i3 

7.4 X loll 

1.5 X 1013 

(3.07,3.37) 

3.22 

1.2 X 

IO13 

1.8 X lOii 

3.3 X IO13 

7.7 X 10i3 

3.9 X 10i3 

1.2 X 10i3 

3.1 X lOii 

7.7 X 10i3 

(2.81,3.07) 

2.93 

8.9 X 

loll 

1.2 X loll 

2.1 X 1013 

4.5 X 1013 

2.2 X 1013 

6.5 X IQii 

1.6 X loll 

4.5 X 1013 

(2.59,2.81) 

2.70 

6.5 X 

loll 

8.2 X IO13 

1.3 X IO13 

2.7 X 10i3 

1.3 X 10i3 

3.6 X IQii 

8.0 X 10i° 

2.7 X 10i3 

(2.44,2.59) 

2.52 

5.0 X 

loll 

5.9 X IO13 

9.2 X 10i3 

1.8 X 10i3 

8.0 X lOii 

2.1 X lOii 


1.8 X 10i3 

(2.28,2.41) 

2.34 

3.7 X 

loll 

4.2 X IO13 

6.0 X 10i3 

1.1 X 10i3 

4.8 X lOii 

1.2 X lOii 


1.1 X 10i3 

(2.19,2.28) 

2.23 

3.0 X 

loll 

3.3 X 1013 

4.6 X I012 


3.4 X loll 

8.2 X 10i° 


8.0 X loll 

(2.11,2.19) 

2.15 

2.6 X 

loll 

2.7 X IO13 

3.6 X 10i3 


2.5 X lOii 

5.9 X 10i° 


6.1 X lOii 

(2.03,2.11) 

2.07 

2.2 X 

loll 

2.2 X IO13 

2.9 X 10i3 


1.9 X IQii 



4.7 X lOii 

(1.96,2.03) 

2.00 

1.9 X 

loll 

1.8 X IO13 

2.3 X 10i3 


1.4 X IQii 



3.6 X lOii 

(1.87,1.92) 

1.89 

1.5 X 

loll 

1.4 X 1013 

1.6 X 1013 


9.2 X IQil 



2.5 X loll 

(1.77,1.81) 

1.79 

1.2 X 

loll 

1.0 X 1013 

1.1 X 1013 


5.8 X IQil 



1.6 X loll 

(1.69,1.72) 

1.70 

9.2 X 

IO13 

7.5 X 10i3 






1.1 X lOii 

(1.53,1.56) 

1.55 

5.8 X 

IO13 

4.2 X 10i3 






4.9 X lOil 

(1.30,1.32) 

1.31 

2.5 X 

IO13 

1.5 X I0I2 







(1.12,1.14) 

1.13 

1.1 X 1013 









Table 1. Column 1: The v{M, z) bins of halos used to calibrate the non-linear bias. Column 2: The corresponding average peak height, 
Vav Other columns give the average mass, Mav, of halos in the sample at the given redshift. If Mav is not given, then the v(M,z) bin 
is not used in our analysis. Columns 3-6 results are for MXXL halos, columns 7-9 are for MS-W7 halos whereas column 10 refers to 
Millennium halos. 


halos in different mass bins. However, it would be more use¬ 
ful if the non-linear bias conld be calcnlated from the dark 
matter power spectrnm. One could then hope to apply our 
results in more general contexts. Furthermore, as discussed 
in the previous section, we first focus on high redshift ha¬ 
los. Therefore, we measured the scale-dependent halo bias 
{hni{r, M, z)) of dark matter halos from the MS-W7, MXXL 
and Millennium simulations in bins of peak height, i^(M, z), 
in the redshift range 2 — 5. 


First, we briefly discuss the effects of resolution and 
the halo exclusion effect on the estimated non-linear bias by 
comparing results obtained using MXXL and Millennium 
halos at 2 = 3. As discussed before, these simulations have 
the same set of cosmological parameters, but different mass 
resolutions and volumes. In Fig.[^ we have plotted the non¬ 
linear and scale-dependent bias at z = 3 for halos (in a given 
iz bin) from the MXXL and Millennium simulations. Also 
shown by the vertical line, is the length scale corresponding 
to twice the virial radius of the most massive halo in each 
sample. 


We first note that in Fig. the scale-dependent bias 
at 2 : = 3 is expressed as a function of z) at the same 

redshift. This choice has been made in several analytic and 
numerical studies probing the clustering of dark matter ha¬ 
los on quasi-linear scales | Hamana et al. pooI| |Scannapieco| 
&_^axkanajj2002| JDiaferio_et_^^ 200^ Tinker et al.||2005[ 


Reed et al. 12009 van den Bosch et al. 112013 1. These studies 


present the scale-dependent bias, bni (r, M, z) as a function of 


the non-linear dark matter correlation function, z) at 

the same redshift. In what follows, we adopt this approach 
and express b„i(r, M, z) as a function of CmnK^i^)- In this 
approach, the scale-dependent bias can be thought of as a 
function of the non-linear dark matter power spectrum, since 
the matter correlation function is the Fourier transform of 
the non-linear dark matter power spectrum. 

It is clear from Fig.[^that the halo bias estimated from 
the two simulations agree well with one another on scales 
larger than twice the virial radius of the most massive halo 
in the sample. However, on smaller scales, the estimated halo 
bias is different between the two simulations. Since both sim¬ 
ulations use the same set of cosmological parameters, this 
could be due to the difference in mass resolution between 
the simulations. We also note that, on the largest scales 
(Cmm (^5 2 ) < 0.05), the bias is approximately a constant. On 
smaller scales, the non-linear bias increases with decreasing 
scale (increasing ^iOm) and reaches a maximum value around 
the scale corresponding to twice the virial radius of the most 
massive halo in the sample. On smaller scales than this, the 
halo bias drops to 0. This suggest that, while probing the 
clustering of a sample of halos, the halo exclusion effect is 
important on scales smaller than twice the virial radius of 
the most massive halo in the sample. Because of these effects 
due to the resolution and halo exclusion, our further discus¬ 
sion and analysis will consider the clustering of a sample of 
halos only on scales larger than twice the virial radius of the 
most massive halo in that sample. 
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We now present our estimates of the non-linear bias 
from different simulations at various redshifts in Fig. as a 
function of 2 ). The results measured from the MXXL 

simulation at z = 2 are shown by solid black lines and at 
z = 3 are shown using red circles. All the other curves at 
2 ; = 3,4 and 5 are estimates of b„i{r, M, z) from the MS- 
W7 simulation. Each panel corresponds to a different bin of 
u{M,z) (see labels). We again emphasize that bni{M,r, z) 
at a given redshift is plotted against at the same 

redshift. 

We first note that, on scales corresponding to < 
0.1, b„i(r, M, z) measured from the MS-W7 and MXXL sim¬ 
ulations agree well with each other. On smaller scales, the 
estimated bias is different for the two simulations. The halo 
bias measured from the MXXL simulation drops to zero on 
larger scales compared to the bias estimated from the MS- 
W7 simulation. This is because, for halos in a given v bin, 
the masses and virial radii of MXXL halos are larger than 
those of MS-W7 halos. 


3.2 A model for the halo bias 

As discussed before, b^i (r, M, z) at different redshifts is fairly 
constant for < 0.05. These scales typically correspond 
to comoving length scales greater than 10 h~^ Mpc. Thus 
on such large scales, the expression for the non-linear bias 
reverts back to the usual scale-independent large scale bias, 
which is only a function of the peak height v alone. There¬ 
fore, one can write 

bni{r,M,z) = ^{r,M,z)b{v), (11) 

where b{i>) is the large scale bias. Here the non-linear bias, 
bni{r, M, z), is written as the product of a scale-dependent 
function, 'yir, M, z), and the large scale bias. The scale- 
dependent function 7 (r, M, z) is thus expected to be close 
to unity on large scales. 

To understand the evolution of the non-linear bias of 
rare halos with redshift, one has to calibrate the expres¬ 
sions for the large scale linear bias b{v) and scale dependent 
function 7 (r, M, z). In what follows, we first obtain a fitting 
function for b{v) and then constrain the functional form of 
7 (r, M, z). 


3.3 The large scale bias 


To estimate b{i'), we measured the correlation functions of 
dark matter halos in different v{M, z) bins in the redshift 
range 0 — 5. These :/-bins are given in column 1 of Table 
Also given in the table are the average peak height, Vav and 
average mass of halos, Mav, in these bins at each redshift. 
The average peak height is given by Vav ~ bcjoav, where 
(Tav is computed as 


a 


2 

av 


(^) 


f dM n{M,z) a^(M,z) 
/ dM n{M, z) ' 


( 12 ) 


We define the large scale bias, b{y), as the average bias 
of halos that are separated by 10 ^ r ^ 25 Mpc/h, i.e. 


b{u) = 


\ 


We then obtain a fitting function for b(v) by refitting the free 



J_I_I_L 

12 3 4 

v{M) 


Figure 4. Our fit to Eq. 0 for the large scale linear halo bias 
(solid black curve) along with the simulation measurements (sym¬ 
bols) at various redshifts. The points for 2 = 0 and 2 are mea¬ 
sured from the MXXL simulation and those at 2 = 3 and 5 are 
obtained from the MS-W7 simulation. The red dotted and blue 
dashed curves are the fitting functions of the linear halo bias given 
by jSheth &: Tormen| | |199^ (ST) and |Tinker et alT] l |2010^ (Tinker) 
respectively. 


parameters of E q. (|4| to the b(y) measured from the simula¬ 
tions using Eq. ( |13[ ). The best fit parameters are estimated 
to be A = 1.0, a = 0.36,5 = -1.156,5 = 2.18, C = -0.749 
and c = 2.18, treating all bias values with equal weight. 

In Fig.[^ our fit (black solid line) for h{v) is overplotted 
with the symbols measured directly from the simulation at 
different redshifts. The data points at 2 = 0,1, and 2 are 
measured from the MXXL simulation and those at 2 = 3, 4, 
and 5 are obtained from MS-W7 simulation. It is clear from 
Fig. I^that our fitting function for b(u) agrees very well with 
the measurements from the simulations. In fact, we find that, 
the overall agreement between the simulation results and the 
fitting function is within 3%. Also shown in Fig. [ 4 ] a re the 
fitting functions for halo bias from|Sheth fc Tormen] ( |1999[ ) 
(ST, red dotted curve) and Tinker et al. ( 2010[ ) (Tinker, 
blue dashed curve). One can see from the hgure that when 
v(M) < 2, our formula compares well with these two fitting 
functions (particularly with the Tinker formula). However, 
for larger values of v{M), the ST formula predicts lower bias 
values and Tinker formula gives slightly larger bias values 
compared with our formula for halo bias. Thus for rarer ha¬ 
los with V ^ 2, our analysis predicts a slightly lower value for 
the large scale bias compared to the Tinker formula. How¬ 
ever, we also note that Tinker et al.|||2010 1 use the spherical 
overdensity (SO) algorithm ( Tinker et al.|[2008 l to identify 
halos, which is different from the FOF(0.2) algorithm used 
in this work. Such a difference in the bias of halos identified 
by these two algorithms has already been noted by 


et al. (2008 

20101. Further, the simulations used by 

Tinker 

et al. (2010 

1 span a wider range of cosmological parame- 


Tinker 
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z 

Figure 5. The effective power law index, am, (Eq. ( |14[ |) is plot¬ 
ted as a function of redshift. The blue triangles and red circles 
correspond to the values obtained for the MS-W7 and the MXXL 
cosmology respectively. 


ters than used in this work. This could also account for the 
difference between the estimated large scale halo bias. 


3.4 The scale dependence of halo bias 

Having obtained the expression for the large scale (r ^ 10 
Mpc/h) bias of rare halos in the redshift range 0 — 5, we now 
wish to calibrate the scale dependence 7 (M, r, z) of b{r, M, z) 
on quasi-linear scales. We first concentrate on the expression 
for 7 (M, r, z) of halos in the redshift range 2 — 5. This is 
because, as we shall see later, 7 (M, r, z) probably has an 
explicit dependence on the effective matter density Qm{z) of 
the universe as a function of redshift. We expect to separate 
out this dependence by focussing on high redshifts {z ^ 2) 
where 0 .m{z) ~ 1. 

It is clear from Fig. that the scale dependent, non¬ 
linear bias of halos (from the MS-W7 simulation) of a given 
v{M,z) as function of and at z = 3,4 and 5, 

agree fairly well with each other. In fact, in this case, the 
agreement between estimates of b(r, M, z) at different red- 
shifts is better than 10 % on quasi-linear scales (r ^ 15 
h~^ Mpc). However, Fig.j^shows that the scale dependence 
b{r, M, z) of MXXL halos for the same v[M, z) at z = 2 
and 3 (which is also expressed as function of Cmm(?’,-z)) is 
quite different from that of MS-W7 halos at higher redshift. 
Thus, 'y{M,r, z), which accounts for the scale dependence 
of b{r,M,z), cannot be described as a function of just two 
variables, ^mm(r, z) and !^(M, z). This suggests that the non¬ 
linear bias on quasi-linear scales is not a simple function of 
the dark matter power spectrum and any htting function 
should be a function of other parameters. 


Such an explicit dependence of the halo bias on param¬ 
eters other than the dark matter power spectrum has been 


discussed in several analytical studies (|Matsubara 

1999 

IBlanton et al.|19991 |Iliev et al.|20031 |Sheth & Tormen 

2004 


Gao et al.|2005al|Jeong fc Komatsu|2009] [McDonald fc Roy| 


2009 Lazeyras et al.|2016 1. These studies point out that one 


may potentially require an infinite number of parameters to 
express the scale-dependent bias. On the other hand, most of 
the available fits to the results of N-body simulations present 
the scale-dependent halo bias as a universal function of the 


Diaferio 


et al.| 

2003| |Cen et al.||2004| |Tinker et al. 2005 Gao et al. 

2005b 

Reed et al. 

to 

0 

0 

_ 

Desjacques et al. 

20101. We will 


now investiagate whether the scale-dependence of the halo 
bias measured from the simulations can be expressed as a 
function of additional parameters along with and 

v{M,z). 

We find that adding one more parameter can account 
for all the simulation results in the redshift range 2 — 5. That 
is, 7 (M, r, 2 ) at 2 ^ 2 ^ 5 can be expressed, to sufficient 
accuracy, as function of three variables, u{M,z), 
and amiz), an effective power law index of a{M,z). This 
effective power law index is defined as 


Ctm{z) 


log(1.686) 

l0g[M„i(2)/Mco;(2)] ’ 


(14) 


where the non-linear mass scale, M„i{z), and the collapse 
mass scale, Mcoi[z), at any redshift are masses at which the 
peak heights are, respectively, 1.686 and 1. The parameter 
am{z) can be thought of as an effective power law index of 
cr(M, 2 ) in the mass range from the collapse to the non-linear 
mass scale (see Appendix A for more details). The depen¬ 
dence of 7 (M, r, 2 ) on am{z) can be tentatively understood 
from Fig. where we have plotted this ratio as a function 
of 2 . The blue triangles at 2 = 3,4 and 5 are obtained for 
the MS-W7 and the red circles at lower redshifts are for the 
MXXL cosmological parameters. The figure clearly shows 
that am{z) is nearly constant in the redshift range 3 — 5 for 
the MS-W7 cosmology. However, at 2 = 2 where the MXXL 
cosmology is used, am{z) is larger. Such a difference is per¬ 
haps related to the departure from the universal nature of 
non-linear bias as a function of v{M, 2 ) and z). Moti¬ 

vated by this, we further investigated whether 7 (M, r, 2 ) can 
be expressed as a function of u{M,z), Cmm(^)-z) and amiz) 
and we hnd that, it is indeed possible to obtain a good ht 
for high-(j halos (iz > 1). The resulting htting function is 
given by 


7(Cmm,'^, am) = 

(^1 -I- Koil + k3/am)yog (^1 -I- X 

-I- I/o(l + ^s/Om) log -I- Cmm ) ) (1^) 

The free parameters are estimated to be Kq = —0.0697, ki — 
1.1682,1=2 = 4.7577,^3 = -0.1561, Lo = 5.1447, H = 
1.4023, h = 0.5823 and I 3 = —0.1030. In htting Eq. (151, we 
have used the non-linear bias estimated from simulations for 
all v'iM, 2 ) bins in redshift range 2 = 2 — 5 given in Table 
For htting 7 ( 5 ™, j/, am), we considered halo correlation 
functions only on scales larger than twice the virial radius 
of the most massive halo in the sample. Moreover, we have 
restricted our analysis to r ^ 30 Mpc. 

It is also important to note from Tablej^that the masses 
of the halos used in our analysis range typically from 5 x 10^° 
to 5 X Mq. Thus, the high-a halos of interest are 

those expected to host galaxies in the r edsh ift range 2 — 5. 

The expression for 7 (^mm, >z) in Eq. (151 is plotted (solid 


© 0000 RAS, MNRAS 000, 000-000 

















































The clustering of dark matter halos: scale-dependent bias on quasi-linear scales 9 



tsim 
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Smm 


Figure 6. The scale dependence of non-linear bias, measured from N-body simulations, in the redshift range z = 2 —5 

as a function of halos in the i^-bins listed in the legend. Solid lines: the fit for , i', ctm) presented in this work. Dotted 

lines: the fitting function for 7 given by [Reed et al.||2009||. 


line) in Fig. at different redshifts as a function of 
at that redshift and for different values of v{M). In Fig. 
the results shown at z = 2 are from the MXXL simulation 
and those at other redshifts are from the MS-W7 simula¬ 
tion. From Fig. it is c lear that the parameterization of 
7 (^“™, iz, (Jeff) using Eq. (151 fits the whole range of data 


measured directly from the simulations very well. In partic¬ 
ular, we note that our fit is consistent with the results from 
simulations to within an overall accuracy of 4%. This sug¬ 
gests that, it is indeed possible to find a fitting function for 
the scale dependence of the non-linear bias in the redshift 
range 2 — 5, through u and am- 

In Fig.[^ we have also plotted in dotted lines, the fitting 
function for 7 given by | Reed et al.| ( |2009[ ). These authors 
parameterized the scale dependence as a function of the large 
scale bias b(y) and a{r, z) as 


7 ( 5 ( 1 /), (T(r, z)) = [1 + 0.036®(i/)(T^(r, z)]. 


(16) 


It is clear from Fig.j^that the Reed et al. (20091 fit compares 
reasonably well with results from the MXXL and MS-W7 
simulations at z = 3,4 and 5, especially at 2 = 5. However, 
their formula is not quite consistent with the MXXL simula¬ 
tion results at 2 = 2. This is expected, since the |Reed et al.| 
(20091 expression for non-linear bias of a halo of mass M at 
a scale r depends only on the dark matter power spectrum 
through the rms linear density fluctuations on the mass scale 
M and length scale r. As noted before, such a simple depen¬ 
dence cannot accurately account for the scale dependence of 
the non-linear bias seen from simulations. 


3.5 The evolution of 7 (M, r, 2 ) to low redshifts 

Having obtained a fitting function for 7 (r, M, 2 ) for high-z 
halos in the entire redshift range 2 — 5, we now include low-z 
data to probe the evolution of the scale-dependence of non- 
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Figure 7. Our fitting function for the scale dependence of the non-linear bias, ™ the redshift range 0 — 5 as a function 

of various choices of v (see legend) is shown using solid lines. The N-body simulation from which results were measured along 

with the redshift is labelled on each panel. 


linear bias from 2 = 0 to 5. We note from Table that, at 
2 = 0 and 1, the masses of the rare dark matter halos used in 
our study ranges from 10^^ to correspondingly 

they host galaxies as well as groups and clusters. 

In Fig. we show 7 (r, M, z) estimated from the simu¬ 
lations (symbols) over the full redshift range. The data at 
2 = 3,4 and 5 are measured from MS-W7 and those at 
2 = 0,1 and 2 are from the MXXL simulation. It is clear 
from the figure that at lower redshifts (2 = 0 and 1) the 
scale dependence of the non-linear bias is rather weaker com¬ 
pared to other redshifts. Such a weak scale-dependence of 
halo bias on quasi-linear scales at lower redshifts (z = 1) 


can be found also in the the analytic work of |Scannapieco| 


& Barkana (20021. In particular, at 2 = 0, the halo bias 
increases by only ~ 10% on quasi-linear scales even for the 
most massive and hence rarest halos at that redshift. 

It turns out that one can obtain a fit for the non-linear 
bias which extends to redshift 0 by adding an additional 
parameter, flm(z), the matter density of the universe at a 
given redshift. 

Ilm (1 + 


^m(zj — 


IIm(l + z)^ -(- Qa 


(17) 


Thus, the evolution of 7 (r, M, z) in the redshift range 0 — 5 
can be expressed as a function of four variables, flm(z), 
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v{M,z), CA'tt{r,z) and a^{z). In particular, we obtained 
a fitting function for 'y{r, M, z) using the non-linear bias 
estimated from simulations for halos in bins of v{M, z) in 
the redshift range z = 0 — 5 given in Table As before, 
we have used the correlation functions only on scales larger 
than twice the virial radius of the biggest halo in the sample 
and smaller than 30 Mpc/h for the analysis. The resulting 
fitting function is given by, 

(^1 -f A:o( 1 -f kslOm) {^m{z))^'^ log (^1 -I- ) X 

+ Lo{l + Iz/am) {0,-m{z)y* log -I- ^ 

(18) 

Here Kq = 0.1699, fci = 1.194, fc 2 = 4.311, fca = 

-0.0348,^4 = 17.8283, Lo = 2.9138, H = 1.3502, ia = 
1 . 9733,^3 = —0.1029 and U = 3.1731. The fitting function 
in Eq. ( |18[ ) is plotted as solid lines in Fig. along with data 
points measured from simulations. Our fit is in remarkable 
agreement with data from all the simulations over the en¬ 
tire range of redshifts from 0 to 5, peak heights and length 
scales. The overall agreement of this fit with the data from 
the simulations is found to be better than 4%. 


3.6 Halo clustering as a function of the linear 
matter correlation function 


We have, so far, presented a model for the non-linear cluster¬ 
ing of dark matter halos as a function of the non-linear dark 
matter correlation function, measured from the 

simulations. In this section, we model halo clustering as a 
function of the linear matter correlation function, z). 

This is well motivated because z) is easier to compute 

without uncertainties, compared to the non-linear matter 
correlation function. Thus, for all practical purposes, it will 
be convenient to express the non-linear bias as a function 
of 5mm (^1 • 2 ). The linear matter correlation function is com¬ 
puted from the linear matter power spectrum as 

00 

eZ.{r,z) = I P‘*-{k,z) (19) 

0 

In order to model the non-linear halo bias as a function 
of 5mm (’'1 z), we first define b„i{r, M, z) at any given scale as 


6 „i(r, M, z) 



( 20 ) 


The new definition of bni (r, M, z) is similar to that given by 
Eq. (10 1, but uses 5mm(*'> 2 ) instead of 5mS)(’'i-2)- Following 
section |3)2| we then express the non-linear bias as the prod¬ 
uct of the scale-independent large scale bias b{u) and the 
scale-dependent function 'y{r, M, z) (see Eq. 0 ). A new 
fitting function is obtained for the large scale bias by re¬ 
fitting the free parameters of Eq. ||^ to the large scale bias 
measured from the simulations. The new best fit parameters 
are given by A = 1.0, a = 0.223,5 = 1.156,6 = 2.167, C = 
-0.748 and c = 2.167. 


As before, we find that, the scale-dependence of halo 
bias 7 (r, M, z) can be expressed as a function of v, 
am,^m{z) and the linear matter correlation function, 5mm- 


The fitting function f or 7 is assumed to have the same func¬ 
tional form as in Eq. (18 1 with 5mm being replaced by 5mm- 
The free parameters of Eq. (18 1 are then determined by 
fitting this equation to the data measured from all the sim¬ 
ulations in the redshift range 0 — 5. The new best fit pa¬ 
rameters of Eq. (18 1 are given by Kq — 0.000529, ki — 
1.0686,^2 = 3.4158,^3 = -204.1715,^4 = 26.9453, Lo = 
0.448,/i = 2.128,^2 = 3 . 0222,13 = 0.226 and h = 1.691. We 
emphasize that the new fit for 7 as function of 5 mm agrees 
very well with the simulation data. The overall agreement of 
this fit with the data from all the simulations given in Table 
[^is found to be better than 5%. 


4 DISCUSSION AND CONCLUSIONS 


We have revisited the problem of modelling the non-linear 
clustering of rare dark matter halos, that collapse from 1 —3 (j 
fluctuations, on quasi-linear scales. In particular, we found 
using high-resolution N-body simulations that the non-linear 
bias of high redshift galactic dark matter halos is strongly 
scale dependent on scales ~ 0.5 — 10 h~^ Mpc. These scales, 
commonly referred to as quasi-linear scales, correspond to 
scales larger than the typical virial radii of dark matter ha¬ 
los. Even though we primarily focussed on the clustering of 
dark matter halos in the redshift range 0 — 5, our results 
are applicable to higher redshifts, including the cosmic dark 
ages before the epoch of reionization. 

First, we estimated the correlation functions of dark 
matter halos at z = 2,3,4 and 5 from the N-body simula¬ 
tions, in mass bins in the mass range 10^^ — 4 x lO^^M©. 
These are the typical masses of dark matter halos that host 
LBGs and LAEs in the same redshift range and correspond 


to rarer objects collapsing from high cr fluctuations (Jose 
et al.|[^13b l. We then showed that, on quasi-linear scales, 
there is a strong discrepancy between the halo correlation 
functions computed using the scale independent, linear halo 
bias and those measured directly from simulations. This sug¬ 
gests that the linear bias approximation is not sufficient to 
explain the clustering of high-z, rarer dark matter halos on 
quasi-linear scales. 

To quantify the non-linear bias of dark matter halos in 
this redshift range, we measured the correlation functions of 
halos, from simulations, in bins of halo peak height, iz(M, z). 
The non-linear bias is defined as the square root of the ratio 
of halo and dark matter correlation functions (see Eq. ( |10[ )). 
We found that the non-linear bias of a halo can be expressed 
as the product of the usual scale independent large scale bias 
6 (M, z) and a scale dependent function 7 (r, M, z) (Eq. 0)- 
We also obtained a fitting function for b{M, z) which de¬ 
pends only on the peak height, !z(M, z), of dark matter ha¬ 
los. This fit compares very well with other formulae for large 


scale bias in the literature (Sheth & Tormen 1999 Tinker 
et al.| 2010 | , especially for halos with v <2, collapsing from 
low (7 fluctuations. For rarer halos with larger values of v, we 
obtained a slightly lower value for the large scale bias com¬ 


pared to the formula given by Tinker et al. (20101. However, 


as noted before, this could be due to the difference between 
the SO and FOF(0.2) halo finder algorithms respectively 


used by Tinker et al. (20101 and in our work. Further, both 
studies use distinct simulations with different cosmological 
parameters for calibrating the bias. 
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We find that, for 2 = 2 — 5, the scale dependence of 
the non-linear bias, 'y{r, M, z), for halos of mass M at any 
length scale r depends on three parameters, the peak height, 
i'{M, z), of halos at mass M, the dark matter correlation 
function (^mm(>',-z)) at that length scale and am, an effec¬ 
tive power law index of a{M) at that redshift. We obtained 
a fitting function that describes the scale dependence of 
7 (r, M, z) as a function of these parameters in the same red- 
shift range. Our fit agrees with the simulation results within 
an accuracy of 4%. 

The scale dependence of non-linear bias at a sc ale r is 
usually parametrized in real space using 2 ) (Tinker 

et al.|2005 1 or the rms linear overdensity in uniformly over¬ 


dense spheres of radius r, j(r, 2 ) (Hamana et al.|2001 Diafe- 


rio et al.|2003 Reed et al.|2009 (. Both (j{r, 2 ) and z) 


can be expressed as functions of the dark matter power spec¬ 
trum. However, we find that the scale dependence of the 
bias, as quantified in terms of 'yir, M, z), is not described 
by such parametrizations, but rather depends on the quan¬ 
tity am{z). But, to compute am{z), one requires only the 
linear dark matter power spectrum. Therefore, it can be ar¬ 
gued that at high redshifts (2 ^ 2 ) the non-linear bias is a 
universal function of the dark matter power spectrum. 

We extended our analysis by probing the non-linear 
bias of low redshift, rarer dark matter halos on quasi-linear 
scales, using MXXL halos at 2 = 0 and 1. Interestingly at 
lower redshifts, especially at 2 ~ 0 , the scale dependence of 
non-linear bias is weaker than at high redshifts and is within 
10-20 % of the large scale bias measured from simulations at 
any scale. We propose a fitting function for the non-linear 
bias as a function of the matter density of the universe at 
a given redshift {Q,m{z)) along with v{M,z), QJ^{r,z) and 
am{z). Remarkably, this fitting function, calibrated using 
the MS-W7 and MXXL simulations, captures the redshift 
evolution of non-linear bias for a wide range of halo masses 
and length scales within an overall accuracy of 4%. 

The dependence of 7 (r, M, z) on Q.m{z) at low redshifts 
breaks the universality of the non-linear bias with respect to 
the linear matter fluctuation field. Thus the observed large 
scale bias of any galaxy population, which depends only on 
the dark matter power spectrum through v(M,z), will not 
uniquely determine the scale dependence of the bias. This 
may provide an opportunity to use the scale-dependence of 
halo bias as a valuable tool to probe cosmology, particularly 
the matter density of the universe. 

We have also extended our analysis by expressing the 
non-linear bias as a function of the linear matter correla¬ 
tion function, z). Here also the non-linear bias is ex¬ 

pressed as the product of the scale independent large scale 
bias, 6 (M, 2 ), and the scale dependent function, 'yir, M, z). 
We first obtained a fitting function for the large scale bias 
as a function of u{M,z). A htting function for 'y{r,M,z) 
is then obtained as a function of linear matter correlation 
function, along with u(M,z), am(z) and Qm(z). 

The new fit for 7 agrees with the data from the simulations 
within an accuracy better than 5%. We emphasize that this 
model parameterizes the clustering of dark matter halos as 
a function of 2 ) instead of the non-linear matter cor¬ 

relation function. Such a model could be quite useful for 
practical purposes as it is easier to compute z) an¬ 

alytically without uncertainties compared to the non-linear 
matter correlation function. 


In general, the halo bias of high redshift, rare dark mat¬ 
ter halos is significantly non-linear and scale dependent on 
quasi-linear scales. On the other hand, at 2 = 0, this scale 
dependence is quite weak and seems to be dependent on 
the matter density of the universe. The non-linear bias is 
expected to have interesting implications on observations of 
the high redshift universe. For example, the halo occupation 
distribution modelling of LBG clustering at high-z (2 is 3) 
usually assumes a linear halo bias | Hamana et al. 2004 2006| 


Hildebrandt et ar]|2009 Lee et al.||2009 Jose et al.||20131^ . 

However, at these redshifts, the typical LBGs collapse from 
2 — 3a fluctuations. Hence, one has to incorporate the non¬ 
linear bias to improve the the clustering predictions of LBGs 
on quasi-linear scales. Thus the non-linear bias could change 
the predicted shape of the two point correlations functions 
of high redshift LBGs and also LAEs, quasars and even the 
redshifted 21 cm signals from the pre-reionization. It would 
be interesting to explore the implications of the non-linear 
and scale-dependent bias in the high-z universe. To do this, 
one may need to incorporate the effects of baryons on the 
clustering of galaxies through the physics of galaxy forma¬ 
tion and also of assembly bias. This is left for future work. 
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APPENDIX A: THE EFFECTIVE POWER LAW 
INDEX OF a{M) 


For any mass scale M, the variance of smoothed d ensity 
contrast a^{M) oc fcifP(feM)crifc®+"=« |Peebles||l980| . Here 
rieft is the effective spectral index, which is ~ —2 on galactic 
scales and —1 on cluster scales. We also have ~ 

Thus we get 


— (3 + rieff) 


a(M) oc M 6 


oc M 


(Al) 


Given the non-linear mass, Mni, and collapse mass, Mcoi, 
corresponds to the mass scales where a{M) is respectively 
1 and 1.686, it is possible to define an effective power law 
index am as 


cr(Mc, 


r{M„ 


= 1.686 


Thus we have 



Otrn 


log(1.686) 

\0g{Mnl/Maol) 


0.2269 


log 


Mnl 

Mcol 


(A2) 


(A3) 
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